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ABSTRACT 

Based on the mean-field approximation of relativistic a-co-p model with the nonlinear self-interaction of a meson 
and density-dependent coupling constants respectively, the nuclear symmetry energy S, its slope L and curvature K as well 
as the proton fraction xp and the pressure P of the matter in the outer core of neutron star are studied. It is found that 
greater S in high-density region and L at saturation density produce greater proton fraction and pressure of neutron star 
outer core matter. And the density threshold of proton super fluidity and the direct URCA process, which are model 
dependent, decrease as the nuclear symmetry energy increases. 
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INTRODUCTION 

S(p B ) is defined as the nuclear symmetry energy (NSE) if the energy per particle of nuclear matter is expanded as 
the following seriesfl, 2]: 

E(p B > a) = e (p b -0)+ S(p B )■• a 1 + S 4 {p b )■• a 4 +.... (1 ) 

It means the impact of isospin asymmetry a=(p n -p p )/p B on energy per particle of nuclear matter and reflects the 
difficulties of nuclear matter from symmetry to asymmetry. Pis=p n +P p with p„, p p and p B to be the densities of neutron, 
proton, and nucleon, respectively. 

NSE has been studied extensively in many experiments [3], among which the more prominent are the 
measurement of nuclear polarizability[4] and nuclear mass[5], parity violation in electro scattering on nuclei[6], the giant 
dipole resonances[7], heavy ion collisions[8-10] and nucleon optical potential[11], direct detection of gravitational 
waves[12], X-ray observations of neutron stars[13, 14] etc. The theoretical works of NSE mainly include the research of 
finite nuclei and infinite nuclear matter [15-25]. The NSE J ~ 32 MeV[3, 5] at saturation density and the binding energy of 
infinite symmetric nuclear matter is roughly equal to 16 MeV with incompressibility coefficient K 0 = 240 + 20 MeV 
here [26]. 

There are two kinds of results of Hartree-Fock approximation of asymmetric nuclear matter according to the 
density-dependence of NSE around saturation density [19, 27-31]. Neutron skin of neutron-rich nucleus is sensitive to NSE 
under saturation density [32-36]. To highlight this dependence, the neutron skin has been widely studied by potential 
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theory and field theory [21, 37-40]. Similarly, the NSE has an important influence on the properties of a neutron star of 
which the outer core matter consisted mainly of neutrons, protons, electrons, and muons can be simulated by homogeneous 
infinite nuclear matter above saturation density [41-43]. The fraction of particles in a neutron star, and therefore, the 
equation of state (EOS) of neutron star matter are affected by NSE [43-50]. EOS is a key input when we study the 
properties of neutron stars, so the radius of a neutron star [51, 52], tidal deformation [53-55], the core-crust transition [56, 
57], and other gravitational properties [58-61] are all related to NSE. And most directly, the fractions of protons x p are 
constrained by NSE [62]. The 1 S 0 proton pairing is predicted to exist, while the neutron superfluidity may occur mainly in 
the coupled 3 P 2 - 3 F 2 two-neutron channel in the quantum fluid in the outer core when the proton density increases to the 
range about Pq/10. The superfluidity of neutrons, protons, hyperons in the S () channel and neutrons in the 3 P 2 channel has 
been studied extensively [63]. However, density range where energy gaps appear in neutron stars is still unclear. On the 
other hand, the direct URCA process which will be restrained by superfluidity of nucleons is allowed at x p ~ l/9[44, 64- 
66]. The competition between superfluidity of nucleons and direct URCA process plays an important role in evolution of 
neutron stars [67-70]. The direct URCA process is allowed in several models with appropriate density-dependence of NSE 
when the density is over maximum density 0.5 fm' 3 , while it is not allowed in some other models for the density needed is 
greater than the density of the most massive Stella core. A detailed discussion in this aspect is shown in Ref. [71]. 

Following these works, the impact of NSE on pressure and proton fraction in the outer core of neutron star will be 
studied in the relativistic a-co-p model and the results under skyrme model will also be reproduced for comparing in this 


paper. 


The basic theories about relativistic a-co-p model and NSE will be introduced in the following section. And the 
numerical calculation results and analysis will be presented in section 3. The summary will be given in the last section. 

Formalism 

The Lagrangian density in the relativistic a-co-p model is 



( 2 ) 


0> po ='d M V v -d v V /t , 


where i// B , i//;, 0 , V / / „and b,, are the field operators of nucleons (neutron and proton), leptons (electron and muon), 
a , co and p meson, and M, m,, m„, m„, and m p their mass, respectively. g a , g aJi and g p are the coupling constants between 
nucleons and a, oo, p meson, respectively. 

By the mean field approximation, we can obtain the expression of the energy density and pressure of infinite 
nuclear matter: 



(3) 



2 


3 ! 


4 ! 
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2 0 3 ! 4 ! 




p B a 


where 0 , V 0 and b 0 are the mean field of a, co and p meson, respectively, and the nucleon density n = c V £ 


3 

/B ’ 


E b — fk 


fB + m B 2 with effective nucleon mass n i B = m B 
The energy per nucleon is defined by 


. , and El = ^jk 2 fA + ml 2 with effective lepton mass m * = m A 


■E{p B ,a) = —. 

Pb 


(5) 


'S(p B )-J + L sym x + — K sym x~ +...; 3 : - 


_ Pb A) 

3/> 0 


(6) 


The NSE can be expanded as following again: 

where the slope parameter / = i£.| and the curvature parameter v =lid . L svm =58.7+28.1 MeV[3,5], - 

s y m a* I s y m 3j 2 h = o 

400</T svm < 100 MeV [21,72]. 


The coupling constants are treated in two ways to reproduce the saturation properties of infinite symmetric 
nuclear matter. The first is called non-linear model (NL) in which a, co, and p are treated as constants and the nonlinear 
self-interaction of a meson c, d is introduced; and the second is called density-dependent coupling model (DDC) in which 
g a , g m and gp are treated as density-dependent and c=d=0, and the density-dependent formulas of the coupling constants g„, 
g fl)i and g p are assumed to be 

1 +a a {x B + c a ) 2 ( 7 ) 

o a o <t0 , , ( \2 ’ 

\ + b a {x B +c a ) 


S CO S CL 


l + a^B+cJ 2 

1 +b m { x B + cJ 


( 8 ) 


* & p S pQ ’ e 




(9) 


where x b =Pb/po- 

Neutrons will decay into protons in the outer core of neutron star if p n >p p +p e , thus the chemical equilibrium 
condition must be considered[73]: 


Pn=Pp+Pe’ 


( 10 ) 


P,=Pu> 


( 11 ) 


where 


Pi 


',=i 


kjj + m i 


i = n, p,e,ft. 
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Finally, the charge neutrality condition must be considered as well, 

P p =P,+Pu- (12 ) 

Numerical Results 

In the NL model, the parameters are fitted to meet the saturation properties of infinite symmetric nuclear matter as 
listed in Table 1. And in the condition of outer core of neutron star, the fraction of protons x p , the pressure P, NSE S with its 
slope L and curvature parameter K are studied, and the numerical results are depicted in Figure 1. The solid, dashed and 
dot-dashed curves correspond to the cases NL1, NL2, NL3 listed in Table 1, respectively. At a fixed x B , especially in high- 
density region, P and x p increase with S increasing. The values of L and K at saturation density (x ri = 1) are 75-100 MeV 
and -125-125 MeV, respectively, which are similar to previous works. We can also see that greater K 0 produces greater 
values of the above quantities. 

The solid circles in three cases represent the density threshold of proton superfluidity, which are *8=1.93, 2.2, 
2.49 in NL1, NL2, NL3, respectively. It decreases with NSE increasing. The star represents the density threshold of the 
direct URCA process, which is x B =3.95 in NL1 and are larger than 4 in NL2 and NL3. It should be noted that K„ in NL1 
and following DDC1 in which star appear has exceeded a reasonable value, so it is difficult for the direct URCA to take 
place in the outer core. 

In the DDC model, the parameters which reproduce the saturation properties of infinite symmetric nuclear matter 
are listed in Table 2, and the numerical results are displayed in Figure 2. The solid, dashed and dot-dashed curves 
correspond to the cases DDC1, DDC2, and DDC3 listed in Table 2, respectively. The NSE and the pressure in high-density 
region are smaller obviously, while x p is greater than the one in Figure 1. The relations between x p , P and S (or K 0 ) at a 
fixed nucleon density are similar to Figure 1. 

Comparing to Figure 1, the density threshold of proton superfluidity in DDC1, DDC2, DDC3 decrease to 1.89, 
2.13 and 2.41 times p 0 , respectively. The density threshold of the direct URCA process in DDC1 decreases to 3.55 po, and 
are still greater than 4 p 0 in DDC2 and DDC3. The values of L and K at saturation density are 50-80 MeV and -130—70 
MeV, respectively. 

For comparing, the results in the models of Ref. [74] are depicted in Figure 3. The solid, dashed and dot-dashed 
curves correspond to the cases MSL0, GSKi , and SkTl. The NSE, pressure, and x p are the smallest comparing to Figure 2 
and Figure 1. The relations between x p , P and S (or K 0 ) at a fixed nucleon density are similar to the previous two cases. The 
proton superfluidity appears in greater density than in Figure 1, and the direct URCA process does in greater density than 
4 p 0 . The values of L and K at saturation density are 50-75 MeV and -125-125 MeV, respectively, which are very close to 
the results in our works. 

Summary 

We have studied the influence of S, L, K on the fractions of proton and pressure of neutron star outer core matter 
by the meanfield approximation of the relativistic a-co-p model. For comparing, we also reproduce the results in model 
MSL0, GSKD and SkTl. NSE and pressure in NL, DDC, MSL0, GSKL, and SkTl decrease in turn. x p is the largest in the 
DDC model and smallest in MSL0, GSKD, and SkTl. L at the saturation density is about 75-100 MeV, 50-80 MeV, and 
50-75 MeV, while K is -125-125 MeV, -130-70 MeV and -125-125 MeV in NL, DDC, MSL0, GSKD, and SkTl. The 
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density threshold of proton superfluidity in SkTl, GSKI , MSLO, NL and DDC model decrease in turn, while the one of 
the direct URCA processes is greater in NL1 than in DDC1 and are greater than 4 p 0 in other models. The values of these 
quantities are model dependent, but the relations among them are similar in these models, that is, the greater S in high- 
density region and L at the saturation density, the greater x p and P, but the smaller the density threshold of proton 
superfluidity and the direct URCA. In addition, we find that the larger K 0 , the greater S in high-density region. And K,, has 
exceeded the reasonable value when the density threshold of the direct URCA is smaller than 4 p 0 . Therefore, proton 
superfluidity may appear in the outer core, and the direct URCA is more likely to take place in the inner core of neutron 
star. 



Figure 1: The Density Dependence of the Fractions of Protons X p , Pressure P, and NSE S as well as its Slope L and 
Curvature Parameter K. The Solid, Dashed and Dot-Dashed Curves Show the Results of NL1, NL2, NL3, 

Respectively 
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Figure 2: Similar to Figure 1 except the Parameters to be given in Table 2 


Table 1: The Parameters of the Models NL1, NL2, NL3 


Model 

g<T« 

§co0 

§d0 

c 

d 

K 0 

NL1 

10.22 

13.0989 

6.60 

4200 

-327 

385.5 

NL2 

9.88 

12.0346 

6.45 

8200 

-715 

277.7 

NL3 

9.31 

10.6574 

6.90 

12500 

-1100 

221.4 


Table 2: The Parameters of the Models DDC1, DDC2 and DDC3 


Model 

go0 

§to0 

§d0 

a„ 

b„ 

Co 

a*. 

b 0) 

Co 


K„ 

DDC1 

8.85 

11.3112 

6.82 

0.99 

1 

-0.10 

1.05 

1 

-0.12 

0.32 

379.9 

DDC2 

8.05 

9.7499 

7.26 

0.95 

1 

-0.12 

1.06 

1 

-0.11 

0.32 

281.5 

DDC3 

7.23 

8.0697 

7.58 

0.92 

1 

-0.11 

1.10 

1 

-0.11 

0.32 

235.9 
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